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Abstract 

We study the generating function of the spin-spin correlation functions in the ground state 
of the anti-ferromagnetic spin-1/2 Heisenberg chain without magnetic field. We have found its 
fundamental functional relations from those for general correlation functions, which originate in 
the quantum Knizhnik-Zamolodchikov equation. Using these relations, we have calculated the 
explicit form of the generating functions up to n = 6. Accordingly we could obtain the spin-spin 
correlator {S^S z - +k ) up to k = 5. 
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Recently there have been rapid developments in the investigation of the exact calculation 
of the correlation functions for the spin-1/2 Heisenberg chains. Especially for the ground 
state without magnetic field, explicit analytical form of several correlation functions have 
been calculated in the thermodynamic limit. In this letter we shall report further results for 
exact calculation of correlation functions for the antiferromagnetic Heisenberg XXX chain, 
whose Hamiltonian is given by 

oo 

TC = [SJ SJ +1 + SjSj +1 + Sj 5*1+1] , (1) 

j=-oo 

where S" 1 = cr°/2 with cr" being the Pauli matrices acting on the j-th site. The Hamilto- 
nian can be diagonalized by Bethe ansatz and many bulk physical quantities have been 
evaluated in the thermodynamic limit On the other hand, the exact calculation of 

the correlation functions is still an ongoing problem. 

First, let us list the known exact results of the two-point spin-spin correlators (SjSj +k ), 
which are physically most important correlation functions : 

(SjSfa) = ^ - iCa(l) = -0.147715726853315 • • • , (2) 
(S*S* +2 ) = l~ - ^Ca(l) + Ca(3) = 0.060679769956435 ■ • • , (3) 
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<^ + s) = 4 - + £C.(3) - ^C(1)C.(3) - ^Ca(3) 
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" -g-Ca(5) + -g-Ca(l)Ca(5) 

= -0.050248627257235 • • • , (4) 
<^S* +4 ) = 1 - fCa(l) + ^C(3) - 72C a (l)Ca(3) - ^C(3) 2 - ™C(5) 

+ ^aDas) - f c(%(5) - <£W 

455 3Q90 
+ ^Ca(7) - ^C.(l)Ca(7) + 280C(3)C(7) 

= 0.034652776982728 ■ ■ • . (5) 

Here Ca( s ) is the alternating zeta function defined by ( a (s) = J^°=i( — l)™ -1 / 77,8 ' which is 
related to the Riemann zeta function £(s) = £ a (s)/(l — 2 1_s ). Note that the alternating 
zeta function is regular at s — 1 and is given by Ca(l) = hi 2. Thus analytical form of the 
spin-spin correlators (SjSj +k ) have been obtained up to k = 4 so far. The nearest-neighbor 
correlator (J2J) was obtained from Hulthen's result of the ground state energy in 1938 
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The second-neighbor correlator (jSJ) was obtained by Takahashi in 1977 4] via the strong 
coupling expansion of the ground state energy for the half-filled Hubbard chain. After that 
it had taken some long time before the explicit form of the third-neighbor correlator was 
obtained by Sakai Sjnroishi, Nishiyama and Takahashi in 2003 g They applied the Boos- 
Korepin method [6, 7] to eva 



developed by Kyoto group 



re multiple integral formula for correlation functions 



uate t 



9j, llCl [ill, ll2|](see also [l3j]). The Boos-Korepin method was 
originally devised to calculate a special correlation function called the emptiness formation 
probability (EFP) which is defined by 

This is the probability to find a ferromagnetic string of length n in the ground state. The 
multiple integral formula for P(n) as well as other correlation functions among the adjacent 
n-spins consists of the n-dimensional integrals. 

By the use of Boos-Korepin method, we can reduce the multiple integrals to one- 
dimensional ones by transforming the integrand to a certain canonical form. In the case 
of the present XXX model, the remaining one- dimensional integrals can be further ex- 
pressed in terms of Ca{s), s G Z>j. Hence, in principle, we can obtain the explicit form of 
arbitrary correlation functions by Boos-Korepin method. However, practically it becomes a 
tremendously hard task to calculate the canonical form as the integral dimension increases. 
Actually P(5) is the only correlation function which was calculated by Boos-Korepin method 
among those for five lattice sites [l^ . 

To overcome such a difficulty, Boos, Korepin and Smirnov invented an alternative method 
to obtain P(n) Q|. They consider the inhomogeneous XXX model, where each site car- 
ries an inhomogeneous parameter Xj. Accordingly the correlation functions depend on the 
parameters Xj, for example, P n (Xi, A n ). Boos, Korepin and Smirnov derived simple func- 
tional relations, which P n (Xi, A n ) should satisfy from the underlying quantum Knizhnik- 
Zamolodchikov (qKZ) equations. Together with a simple ansatz for the final form of the 
P n (Xi, A n ), they have shown the problem reduces to solving large linear systems of equa- 
tions for rational coefficients. In this way they could calculate P n (Xi, A„) up to n = 6. 
Especially, by taking the homogeneous limit Xj — > 0, they obtained the analytical form of 
P(n) up to n = 6. Thus this new method is more powerful than the original Boos-Korepin 
method to evaluate the multiple integrals. It was recently generalized to calculate arbitrary 



inhomogeneous correlation functions by Boos, Shiroishi and Takahashi [la ]. They have cal- 
culated all the independent correlation functions among five lattice sites and especially the 
fourth- neighbor correlator (jSJ). In this letter we explore these functional methods further. 
Especially main new aspect is to introduce the generating function of the spin-spin corre- 



lators 
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18, 



19 



20 



21] into the scheme, which allows us to calculate the (S , |S , | +fc ) more 



efficiently 

We remark that Boos, Jimbo, Miwa, Smirnov and Takeyama have obtained more explicit 
recursion relations for arbitrary correlation functions (see also 2) by investigating the 
qKZ equations more profoundly. Their formulas have proven the ansatz for the final form 
of the inhomogeneous correlation functions. However, the formulas contain unusual trace 
functions with a fractional dimension, effective evaluation of which has not been developed 
yet. 



Here we define the generating function by 



2 + ^ 



K, 2"^ 



(Ai, A 2 , A n ), 



(7) 



where k is a parameter. Once we obtain the generating functions, we can calculate the 
two-point spin-spin correlators through the relation j^J 

1 d 2 



K S(S%. +1 ) (Ai, Afc+i) 



-Pfc+i(Ai, A 2 , Afc+i) — P£ (Ai, Aa,.) 



2 8k 2 

- Pk(^2, Afc+i) + P^„ 1 (A 2 , A fc ) j 



K = l 



1 

4' 



(S*S* k+1 ) = lim (SfS z k+1 ) (Ax, A fe+1 ). 
Note that -P^(Ai, A 2 , A„) is a natural generalization of the EFP as 



JK=0 



(Ai, A 2 , A n ) — P n (Ai, A 2 , A n ), 



(8) 



(9) 



with other relations 



Pn (Ab ^2, A n ) — 1, 

n 

Pn~ ^2, A n ) = 2 n (J^ Sj)(Xi, A 2 , A r 



From the qKZ equations for the correlation functions 
function satisfies the following functional relations : 
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(10) 
(11) 



231 ] . we find the generating 



• Translational invariance 



P n K (Ax +X,...,\n + X)= F*(Ai, . . . , A n ) (12) 

Negating relation 

P£(\i, • • • , A n ) = P£(— Ai, . . . , — A n ) (13) 

Symmetry relation 

P„{Xi, \j, Xj+i, A n ) = -P^(Ai, . . . , Xj, A n ) (14) 

First recurrent relation 

Pni^li ■ ■ ■ i ^n-lj A n -1 ± 1) = «P*_ 2 (Ai, . . . , A n _ 2 ) (15) 

Second recurrent relation 



^ lim P n K (A 1; . . . , A n _x, A n ) = ii^P^^Ax, . . . , A n _ x ) (16) 



Here and hereafter we follow the notations in |22j for the spectral parameter Xj and also 
the transcendental function lu(X) below. One easily sees, these functional relations reduce 

n 

to those for the EFP, P n {X u . . . , X n ) [1J if we set k = 0. According to the (proved) general 
form of the correlation functions and also by the symmetry relation ()14j) . we can assume the 
generating function in the form 



P n t (Ai, . . . , A 

ifj r i 

^2 S ^,;(Ai, A 2 | ■ ■ ■ \X 2 i-x, A a ||A a +i, • • • , A n ) JJ^(A 2i _i - X 2j ) + permutations \ , (17) 

1=0 { j=l J 



where 



k=l 
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is a generating function of the alternating zeta values 

w(A) = 2(A 2 - 1) X " k U2k + 1) + - . (19) 

fc=0 

Moreover the function o;(A) has the following properties: 

u(ioo) = 0, w(A±l)=a(A)+7 ± (A)w(A). (20) 

where 



3 1 , . A(A i 2) 



Note also that A^(Ai, A 2 | . . . |A 2 i_i, A 2 i||A 2 ;+i, . . . , A n ) are rational functions depending on 
the parameter k with known denominators: 

A-n,iv^ii A 2 | . . . |A 2 /-i, A 2 ;| |A 2 / + i, . . . , A n ) 

_Hk=l^k-l,2kY\ 2 l+l<k<j<n' K kj , . 

— T=f r Vn,|^li A2| • • • |A 2 ^_i, A 2 ;||A 2 £ + i, . . . , A n J, (^ZZJ 

where A fc j = A fc - Aj and Q^(Ai, A 2 | . . . |A 2 /_i, A 2i ||A 2i+ i, . . . , A n ) are some polynomials. 

Our main proposal is to determine A* l (X 1 , A 2 | . . . | A 2 z— 1 , A^| |A 2 / + i, . . . , A n ) and therefore 
the explicit form of the generating functions from the functional relations (fT2^) - (|TBj) together 
with the ansatz (|T7|) . Simple calculations yield that the second recurrent relation (|T6*j) is 
equivalent to the following recursion equation for A^(. 

1 4- K 

Hm A^(Ai, A 2 | . . . |A 2 ;_i, A 2 ^| |A 2 ; + i, . . . , A n ) = — - — A^_ x ; (Ai, A 2 | . . . |A 2 /_i, A 2 /| |A 2 ; + i, . . . , A n _i). 

(23) 

Note that it is easy to see 

AU\\X 1 ,...,X n ) = (^Y . (24) 



2 

for any n. The first recurrent relation ()15|) can not be reduced to a general formula for A* v 
Instead, we write down the equations for A* t corresponding to the first recurrent relation 
up to n = 3: 

1 + K ^ - 3 -A1 1 (X 1 ,X 1 ±1\\) = k, (25) 



(^ir) 3 + "(AiaKiCAi, A 2 ± 1||A 2 ) - ^(A 2 , A 2 ± l||Ai) = k 



(26) 



A^(Ai, A 2 | I A 2 ± 1) + Tt(Ai 2 )A? 1 (Ai, A 2 ± 1| |A 2 ) = 0. (27) 
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Solving these equations, we have calculated the explicit forms of l up to n = 5, which are 
given by 

Q£i(Ai,A 2 ||) = ~(1-k) 2 , Q^ 1 (A 1 ,A 2 ||A 3 ) = ^-^(1 + A 13 A 23 )Q^ 1 (A 1 ,A 2 ||), 
QUK A a ||As, A 4 ) = ^(1 + A 14 A 24 )g^ 1 (A 1 , A 2 | |A 3 ) - ^p(A? 2 - 4), 
Ql 2 {\ x , A 2 |A 3 , A 4 | |) = 1 (l + |) (1 + 3*)(1 - k) 2 + (a 2 2 - (a 2 4 - 

+ (1 + K) H 1 ~ /t)2 (l + A 13 A 24 )(1 + A 14 A 23 ) - K(1 ~ K)2 (1 - A 13 A 24 )(1 - A 14 A 23 ), 
do 18 

1 + K 

Q5,i(Ai, A 2 | |A 3 , A 4 , A 5 ) = — - — (1 + Ai5A 2 5)Q2,i(Ai, A 2 | |A 3 , A 4 ) 

- 1 ^ (1 ~q )4 (A 2 2 - 4) (4+ A 13 A 23 + A 14 A 24 ), 

Q£ 2 (Ai, A 2 |A 3 , A 4 ||A 5 ) = ^(1 + A 15 A 25 )(1 + A 35 A 45 )g« 2 (Ai, A 2 |A 3 , A 4 | |) 

+ ^^- (Ai2 - 4)(A 2 4 - 4) {(1 + k)(1 - k) 2 (A 14 A 23 + A 13 A 24 ) + 5(1 + k + k 2 + k 3 )} 

— — {lO — 2(A 2 2 + A 2 4 + Ai 4 A 23 Ai 3 A 24 ) + (Ai 4 A 23 + Ai 3 A 24 + 2)(AisA 2 5 + A 3 sA 4 5 — 1)} 

do 

(28) 

One can confirm, the previous known results of two-point spin-spin correlators (0),...,© are 
reproduced through the relation (JHJ). In a similar way, we have calculated -P 6 K (Ai, . . . , A 6 ). 
Since its expression is too complicated, we present here only the final two new results of 
correlation functions among six lattice sites. The first one is the two-point fifth-neighbor 
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spin-spin correlator, 



(SjSj +5 ) 



1 25 >/N 800 >/N 1192 >/N>/N 15368 , . , 9 w 4228 w , 

= 12 " Y Ca(1) + ~ Ca(3) " ^-Ca(l)Ca(3) - ^-Ca(3) 2 - 608Ca(3) 3 - — Ca(5) 

+ ^Ca(l)Ca(5) - ^Ca(3)Ca(5) + 3648Ca ( 1 ) Ca (3) ^ (5) - ^(3)^(5) - ^Ca^f 

+ 5f5a(i)«5)" + ^a3)c(5) 2 + ^c( 5) 3 + g -f U 7) - ^aDC(7) 

11648D 3^39? 
+ 56952^(3)^(7) - — ^— Ca(l)C(3)C(7) - ^-Ca(3) 2 Ca(7) + 7840C Q (5)Ca(7) 

66640 

- 8960Ca(3)Ca(5)Ca(7) - — g— Ca(7) 2 + 31360Ca(l)Ca(7) 2 - 686Ca(9) + 18368C a (l)Ca(9) 

- 53312Ca(3)Ca(9) + 35392C a (l)Ca(3)Ca(9) + 16128C„(3) 2 Ca(9) 
+ 38080Ca(5)Ca(9) - 53760C a (l)Ca(5)Ca(9) 

= -0.030890366647609 ■ ■ • . (29) 
Another one is a six-spin correlation function 

6 

lim Pr^(Ai,...,A 6 )=2 6 (Jj5?> 

= i - 12U1 ) + If CO) - ^(1)C.(3) _ 12^(3)* - ^? Ca (3)3 - ™?<.(5) 
+ !™*<.(1)<.(5) + ^(3)«5) + ^aDU3)U5) - U3)^(5) 
_ ™ Ca(5)2 + "^.(DUS). + 15£* U 3)<.(5)' + i=C»(5)^ + U 7) 
_ ^ C .(l)«7) + 241888C43)C(7) - ^C.(l)«3)«7) - ^«3)'«7) 
+ H^i Co(5)C(7) _ ^ C „ (3)U 5, U 7) - ^«7) 2 + ±^«D«7)* 

_ ^ c „ (9) + ^ U1)m _ ^U3)c (9 ) + ^c.(Da(3)a(9) 

+ »„(3)*f„(9) + ^(5)C(9) - ^ U DC„(5)C(9) 
5 3 3 

= -0.440301669702626 • • • . (30) 

One can see that these analytical expression have the same structure as P(6) previously 
obtained in [l^. Namely they are the polynomials of Co(l) 3 • • • j Co(9) of degree 3 with 
different rational coefficients. We have confirmed these results (J29|) and (|30|) by comparing 
with the numerical data by the exact diagonalization for the finite systems. 
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In summary we have presented a new effective method to calculate the analytical form 
of the spin-spin correlators starting from the generating function. In particular we could 
obtain the fifth-neighbor correlator (SjSj +5 ). The details of the calculations and the further 
results will be reported in a separate publication. We are grateful to K. Sakai, M. Taka- 
hashi and Z. Tsuboi for valuable discussions. This work is in part supported by Grant-in- 
Aid for the Scientific Research (B) No. 14340099 from the Ministry of Education, Culture, 
Sports, Science and Technology, Japan. MS is supported by Grant-in-Aid for young scien- 
tists No. 14740228. 
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